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Abstract. The scaling properties of post-mortem fracture surfaces of brittle (silica glass), duc- 
tile (aluminum alloy) and quasi-brittle (mortar and wood) materials have been investigated. These 
surfaces, studied far from the initiation, were shown to be self-affine. However, the Hurst exponent 
measured along the crack direction is found to be different from the one measured along the propaga- 
tion direction. More generally, a complete description of the scaling properties of these surfaces call for 
the use of the 2D height-height correlation function that involves three exponents £ ~ 0.75, (3 ~ 0.6 
and z ~ 1.25 independent of the material considered as well as of the crack growth velocity. These 
exponents are shown to correspond to the roughness, growth and dynamic exponents respectively , 
as introduced in interface growth models. They are conjectured to be universal. 
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1. Introduction 

Since the early work of Mandelbrot et al. (1984), the study of the morphology of 
fracture surfaces is a very active field of research. Crack surfaces were observed to 
display some universal scaling features although they result from a broad variety of 
specific processes taking place at the microstructure scale (see Bouchaud, 1997, for a 
review). Many experimental results suggested indeed that fracture surfaces are self- 
affine over a wide range of length scales In other words, the height-height correlation 
function Ah(Ar) =< (h(r + Ar) — h(r)) 2 >^ 2 computed along a given direction is 
found to scale as : 

Ah ( Ar\ H 



I v I 

where H is the Hurst exponent and I the topothesy or the scale at which Ah is equal 
to Ar. This exponent was found to be weakly dependent on the nature of the material 
and on the fracture mode with H ~ 0.8 (Bouchaud et al., 1990; Maloy et al., 1992; 
Schmittbuhl et al., 1995; Daguier et al., 1997). This exponent was then conjectured 
to be universal (Bouchaud et al., 1990; MalOy et al., 1992). 
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Since the early 90s, a large amount of theoretical studies suggested scenarios 
to explain these experimental observations. They can be classified into two main 
approaches: 

(i) Percolation-based models where the fracture propagation is assumed to result 
from a damage coalescence process (Roux and Herrmann, 1990; Hansen and 
Schmittbuhl, 2003). These models lead to isotropic fracture surfaces. 

(ii) Elastic string models that modelize the crack front as an elastic line propagating in 
a random medium (Bouchaud et al., 1993; Ramanathan et al., 1997; Ramanathan 
and Fisher, 1997). The fracture surface corresponds then to the trace left by this 
crack front. As emphasized by Ponson et al. (2006), such approaches lead to 
anisotropic fracture surfaces. 

All these models lead to self affine fracture surfaces, but with varying critical expo- 
nents. The relationship between these predicted exponents and the roughness exponent 
measured experimentally remains controversial. 

We present here an extensive experimental investigation of the scaling properties 
of fracture surfaces in heterogeneous materials. Four different materials with four 
different characteristic length scales of the microstructure were investigated: silica 
glass, aluminum alloy, mortar and wood. They were broken using four different types 
of fracture tests and the fracture surfaces were scanned by three different techniques. 
Fracture surfaces observed in all these different materials/failure modes are shown 
to be self-affine, in agreement with results reported in the literature. However, their 
scaling properties are not isotropic as usually believed but require the use a two- 
dimensional (2D) height-height correlation function for a complete description. 

This 2D description involves two independent scaling exponents, which are the 
Hurst exponents measured along the crack front direction and along the crack prop- 
agation direction. Those will be shown to correspond respectively to the roughness 
exponent ( and the growth exponent (3 introduced to describe the depinning of elas- 
tic manifolds in a random medium (Kardar, 1998; Barabasi and Stanley, 1995). To 
describe fully the scaling properties of the experimental fracture surfaces, a third 
exponent - the dynamic exponent z - must be introduced. As expected for Family- 
Vicsek scaling (Family and Vicsek, 1991), this last exponent is shown to be equal to 
the ratio of the two others, z = (/(3. 

These three exponents are found to vary little with the nature of the material and 
the crack propagation velocity. They are conjectured to be universal. Implications of 
these scaling properties will finally be discussed. 

2. Description of fracture tests and surface scanning 

Silica glass and aluminum alloy are archetypes of brittle and ductile materials re- 
spectively while mortar and wood are good examples of isotropic and anisotropic 
quasi-brittle materials respectively. 

(a) Silica fracture surfaces were obtained by applying a DCDC (Double Cleavage 
Drilled Compression) to parallelepipedic (5 x 5 x 25 mm 3 ) samples under stress 
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corrosion in mode I (see Prades et al. (2004) for details). After a transient dynamic 
regime, the crack propagates at a slow velocity through the specimen under stress 
corrosion. This velocity was measured by imaging in real time the crack tip 
propagation through Atomic Force Microscopy (AFM). In the stress corrosion 
regime, the crack growth velocity can be controlled by adjusting the compressive 
load applied to the DCDC specimen (Bonamy et al., 2006). The protocol is then 
the following: (i) a large load is applied to reach a high velocity; (ii) the load is 
decreased to a value lower than the prescribed one; (iii) the load is increased again 
up to the value that corresponds to the prescribed velocity and kept constant. 
This procedure allows to observe on the post-mortem fracture surfaces constant 
velocity zones clearly separated by visible arrest marks. The topography of these 
fracture surfaces is then measured through AFM with in-plane and out-of-plane 
resolutions of the order of 5 nm and 0.1 nm respectively. To ensure that there is 
no bias due to the scanning direction of the AFM tip, each image is scanned in 
two perpendicular directions and the analyses presented hereafter are performed 
on the two sets of images. These images represent a square field of 1 x 1 (1024 
by 1024 pixels). 

(b) Fracture surfaces of the commercial 7475 aluminum alloy were obtained from CT 
(compact tension) specimens, first precracked in fatigue and then broken under 
uniaxial mode I tension. In the tensile zone, the fracture surface has been observed 
with a scanning electron microscope at two tilt angles. A high resolution elevation 
map has been produced from the stereo pair using the cross-correlation based 
surface reconstruction technique described in Amman and Bouchaud (1998). The 
reconstructed image of the topography represents a rectangular field of 565 x 
405 /im (512 by 512 pixels). The in-plane and out-of-plane resolutions are of the 
order of 2 — 3 /um. 

(c) Mortar fracture surfaces were obtained by applying four points bending to a 
notched beam leading to a mode I failure. The displacement is controlled during 
the test. The length of the beam is 1400 mm and its height and thickness are both 
equal to 140 mm. The topography of the fracture surfaces has been recorded using 
an optical profilometer. The maps include 500 profiles of 4096 points each with 
the first profile close to the initial notch. The sampling step along profiles is 
20 fim. Two successive profiles are separated by 50 fj,m along the direction of 
crack propagation. The lateral and vertical accuracy are of the order of 5 fim. 
A transient regime is also observed by a post-mortem analysis of the fracture 
surfaces. On the first 10 mm of the crack propagation, corresponding to the first 
200 profiles, the roughness of the profiles increases with the distance to the initial 
straight notch. A full description of the roughness in this region of the fracture 
surface is given in Mourot et al. (2005). The present study focuses on the geometry 
of the surface far from the initial notch and the first 200 profiles are therefore 
systematically removed from the maps. 

(d) Fractured wood surfaces were obtained from modified Tapered Double Cantilever 
Beam specimens (TDCB) subject to uniaxial tension with a constant opening rate 
leading to mode I failure (see Morel et al. (2003) for details). The wood species 
used in the study is a Spruce (Picea excelsa W.) which is strongly anisotropic. 
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The crack propagated along the longitudinal direction of the wood. As a result, 
the characteristic length scale of the elementary feature of the fracture surface 
is anisotropic : it is respectively of the order of a mm and of a few tenths of 
micrometer in the longitudinal and transverse directions. These values correspond 
respectively to the length and the diameter of the wood cells. As a consequence, 
the height of the surface has been scanned by an optical profilometer over a 
50 x 50 mm area with a higher resolution in the transverse direction (25 /im) 
than in the longitudinal one (2.5 mm) : this map includes 50 profiles parallel to 
the crack front with 2048 points each. As for the mortar fracture surfaces, the 
maps of the wood fracture surfaces correspond to the zone far from the initial 
straight notch where the roughness is statistically stationary, i. e. approximately 
50 mm from initiation. 



Silica glass Aluminium alloy 




Figure 1. Topographic images of fracture surfaces of silica glass, aluminum alloy, mortar and wood. 
Square fracture surfaces were represented here for the sake of clarity. The a>axis and z-axis correspond 
to the crack propagation direction and the crack front direction respectively. 



In all four cases, the reference frame (e x ,e y ,e z ) is chosen so that e x and e z are 
respectively parallel to the direction of crack propagation and to the crack front. 
Figure 1 shows typical snapshots of the fracture surfaces as observed in silica glass, 
aluminum alloy, mortar and wood. These surfaces display striking visual differences 
: the in-plane (along x or z) and out-of-plane (along h) length-scales of the features 
observed strongly depend indeed on the considered material. They are respectively of 
the order of 50 nm and 1 nm for the silica glass surface, approximately 100 [im and 
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30 for aluminum, and 5 and 0.6 mm for mortar. The fracture surface in wood is 
highly anisotropic: The in-plane sizes of the patterns are respectively 50 and 1 mm 
along the longitudinal (x-axis) and transverse (z-axis) directions and 200 fim for out 
of plane sizes. 

Despite their apparent differences, these surfaces share common scaling properties 
to be discussed in the next section. 



3. Fracture surface analysis: Anisotropic scaling properties 

In order to investigate the anisotropy of the experimental fracture surfaces, the ID 
height-height correlation functions Ah(Az) =< (h(z + Az,x) — h(z,x)) 2 > 1//2 along 
the z direction, and Ah(Ax) =< (h(z,x + Ax) — h(z,x)) 2 > x l 2 along the x direction 
were computed for each material. They are represented in Figure 2 for a fractured 
aluminum alloy sample. 

These curves indicate a clear dependence on the measurement direction although 
all profiles are self affine. This anisotropy is reflected not only in the correlation lengths 
and in the amplitudes, but also in the values of the Hurst exponent: Along the crack 
front, it is found to be £ = 0.75 ±0.03, i. e. fairly consistent with the "universal" value 
of the roughness exponent £ ~ 0.8 reported in the literature (Bouchaud et al., 1990; 
Maloy et al., 1992; Schmittbuhl et al., 1993; Daguier et al., 1997). Parallel to the crack 
front, the measured Hurst exponent is significantly smaller, with (5 = 0.58 ± 0.03. 
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Figure 2. ID height-height correlation functions measured parallel to the crack propagation direction 
and to the crack front for an aluminum alloy fracture surface. The straight lines are power law fits. 
The scaling exponents, corresponding to the slope of these lines, are respectively equal to 0.75 and 
0.58 parallel to the crack front and to the crack propagation. 

The observation of two different scaling behaviors in two different directions of the 
fractured aluminum alloy surface suggests a new approach based on the analysis of 
the 2D height-height correlation function defined as: 
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Figure 3. (a): 2D representation of the height-height correlation function for an aluminum fracture 
surface {Ah is normalized by Ax 13 with (3 = 0.58). (b): Domains of different types of variation of the 
height-height correlation function in the (Az, Ax) plane. The grey zone corresponds to a height-height 
correlation function varying as Ah oc Ax 13 



Ah(Az, Ax) =< [h(z + Az,x + Ax) - h(x, z)} 2 >\ /2 (2) 

This function contains informations on the scaling properties of a surface in all di- 
rections. Figure 3a gives a color scale representation of Ah in the (Az, Ax) plane 
for the aluminum alloy. The function Ah was normalized by Ax 13 and logarithmic 
scales were used for all the axis to emphasize the anisotropy of the power-law scalings. 
This representation clearly demonstrates two distinct behaviors of the 2D correlation 
function depending on the orientation of the vector AB of coordinates (Az, Ax). 

If AB lies inside the grey region in Figure 3b (corresponding to the blue domain 
in Figure 3a), the 2D correlation function scales as Ax@ and does not depend on Az. 
The straight boundaries of this domain in these logarithmic scales indicate that its 
width £ (Figure 3b) increases following a power law £ oc Ax 1 / 2 where z — 1.26. 

In other words, from any given point A of the fracture surface, a domain where 
the 2D correlation function scales as Ax 13 develops over a width Az = £ increasing as 
Ax 1 / 2 (the crack propagates parallel to x). Outside of this domain, the 2D correlation 
function depends only on Az. If Ax =0, the correlation function corresponds to that 
presented in Figure 2 and scales as Ah oc Az^ with £ = 0.75. The latter variation is 
not visible in Figure 3a because of the divergence of the normalization term Ax^. 

These scaling behaviors of the 2D correlation function of the aluminum fracture 
surface can be summed up as follows: 

Ah(Az, Ax) = Ax^f(Az/Ax 1 / z ) 

N fl ifu«c ^ 
where j{u)~< r - r ^ 

J v ' \ vS if u» c 

where the constant c is related to the topothesies l x and l z defined along xand z 
respectively. 

The variations of the correlation functions Ah& x are plotted as a function of Az in 
the insets of Figure 4 respectively for silica glass, aluminum alloy, mortar and wood 
fracture surfaces. This Figure demontrates that the scaling properties described by 
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Figure 4- Normalized 2D height- height correlation function variations with Az for various values of 
Ax for silica glass, aluminum alloy, mortar and wood. The data collapse was obtained from Eq. 3 
using exponents reported in Tab. I. 



Equation 3 apply for all fracture surfaces investigated experimentally. For adequate 
values of [3 and z, it can be seen in the main graphs of that same figure that a 
very good collapse of the curves can be obtained by normalizing the x-axis by Ax 1 / 2 
and the y-axis by Aar. As expected from Equation 3, the resulting master curve is 
characterized by a plateau region and followed by a power law variation of exponent 
(. The exponents (3 and z which optimize the collapse, and the £ exponent determined 
by fitting the large scales regime exhibited by the master curves are listed in Table 
I. The three exponents are found to be C — 0.75, (3 ~ 0.6 and z ~ 1.25, independent 
of the material and of the crack growth velocity over the whole range from ultra-slow 
stress corrosion fracture (picometers per second) to rapid failure (some meters per 
second). They can therefore be conjectured to be universal. 
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The ratio of £ to (3 is given in the fourth column. It is worth noting that the 
exponent z verifies the relation z = (3/C, as expected for Family- Vicsek scalings (Family 
and Vicsek, 1991) classically observed in interface growth processes. 

Table I. Scaling exponents measured on fracture surfaces of silica glass, metallic alloy, mortar and 
wood. f3, z and (/ ' (3 are respectively the roughness exponent, the growth exponent and the dynamic 
exponent z while the fourth column contains the ratio of ( to (3. Error bars correspond to a confidence 
interval of 95 %. 





c 


P 


z 


UP 


silica glass 


0.77 ± 0.03 


0.61 ± 0.04 


1.30 ± 0.15 


1.26 


aluminum alloy 


0.75 ± 0.03 


0.58 ± 0.03 


1.26 ± 0.07 


1.29 


mortar 


0.71 ± 0.06 


0.59 ± 0.06 


1.18 ± 0.15 


1.20 


wood 


0.79 ± 0.05 


0.58 ± 0.05 


1.25 ± 0.15 


1.36 



4. Discussion 

Equation 3 suggests analogies with roughening in interface growth processes. In Ref. 
Ponson et al. (2006), surfaces were generated numerically by solving the Edwards- 
Wilkinson (Edwards and Wilkinson, 1982) and the Kardar-Parisi-Zhang (Kardar et al., 
1986) equations, which describe the dynamics of an elastic line perturbed by a thermal 
noise. In such systems, the roughness developed from a flat profile in the transient 
initial regime is known to be characterized by a ID height-height cross-correlation 
function Ah(Az, t) =< (h(z + Az, t) — h(z, t)) 2 >]J 2 with a time dependence verifying 
(Barabasi and Stanley, 1995): 

Ah(Az,t) = t^g(Az/t 1 / z *) 

, , s ( u c ifu<d ^ 

where q(u) ~ < ., , 
yv ; \ 1 if u> d 

where d is constant and (*, (3* and z* refer to the roughness, growth and dynamic 
exponents respectively, as commonly defined in the framework of elastic manifolds 
propagating in a random environment. These three exponents are not independent 
and are related by z* = (*/(3* (Family and Vicsek, 1991). 

In Ref.Ponson et al. (2006), the 2D height-height correlation function Ah(Az, At) 
is computed for these simulated surfaces in the steady state regime reached at long 
times where the roughness becomes time invariant. It is shown that Ah(Az, At) verifies 
the scaling Equation 3 where the exponents (, (3 and z coincide with £*, (3* and z* 
respectively. In other words, the "universal" exponents £ ~ 0.75, (3 ~ 0.6 and (3 ~ 
1.25 measured on experimental fracture surfaces in the previous section correspond 
respectively to the roughness, growth and dynamic exponents. 

Note that these results differ from those of 

previous measurements in the initial transient regime which reported material- 
dependent dynamic exponents (Schmittbuhl et al., 1994; Lopez and Schmittbuhl, 
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1998; Morel et al., 2003; Mourot et al., 2005). These two last studies focused on 
the fracture surface morphology, respectively of mortar and wood, in the roughening 
regime starting from a straight notch. The surfaces were found to be anisotropic and 
to display an anomalous scaling requiring to introduce a fourth, global roughness, 
exponent. 

Moreover, the dynamic exponents were found to be different from those reported 
in the present study on the same materials and to vary from material to material. In 
other words, there is apparently no relationship between the exponents measured in 
the transient and steady state propagation regimes. 

This observation is quite surprising in view of the predictions of theoretical models. 
In models based on elastic lines propagating in random environments, the physical 
mechanism which rules the time evolution of the line is not changing between the 
transient and the long-time regime. Hence, the exponents are found the same in the 
two regimes (Ponson et al., 2006). In experiments, we would thus expect to get also 
the same exponents in transient and stationary regimes. However, in quasi-brittle 
materials, the physical mechanisms that rule the roughening of a fracture surface devel- 
oping from an initial straight notch (microcracks development, material toughness,...) 
is evolving with time during the transient regime. This may explain the apparent 
discrepancy in measured exponents. 



5. Conclusion 

The scaling properties of fracture surfaces of materials with different characteristic 
length scales of their microstructure have been investigated. Measurements were real- 
ized on brittle (amorphous silica), ductile (aluminum alloy) and quasi-brittle (mortar 
and wood) materials. The surfaces were studied far enough from the initiation of the 
crack so that the roughness distribution could be considered as stationary. It was 
first observed on aluminum alloy surfaces that profiles parallel and perpendicular to 
the direction of crack propagation were both self-affine, but corresponded to different 
Hurst exponents. In order to characterize such surfaces, ID measurements such as 
roughness profiles are too restrictive, and the 2D height-height correlation function 
has been shown to be more suitable. For all types of materials investigated, the 
2D correlation function displays similar scaling properties involving three exponents 
£ ~ 0.75, f3 ~ 0.6 and z ~ 1.25 independent of the considered material as well as the 
crack growth velocity. They are shown to correspond respectively to the roughness, 
growth and dynamic exponents introduced in interface growth models. Moreover, the 
three exponents are not independent: z = (/(3 as expected for Family- Viseck scaling 
(Family and Vicsek, 1991). 

The fact that z ^ 1 and (3 ^ C shows that scaling features of fracture surfaces result 
from dynamic processes. In others words, static models like percolation-based mod- 
els (Hansen and Schmittbuhl, 2003) are not suitable to describe fracture surfaces in 
heterogeneous materials. The experimental observations reported here provide strong 
arguments in favor of models based on elastic manifolds moving in a random medium. 
As shown by Roux and Vandembroucq (2002) and Rosso (2002) in simulations of 
the propagation of an elastic line in a random environment, the criticality of such a 
phenomenon is expected to arise only at long times. This is in perfect agreement with 
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the experimental evidence of the existence of three universal exponents in the long time 
regime reported here. It should be emphasized that elastic strings models available 
presently don't succeed in reproducing quantitatively the exponents observed. Work 
in this direction is currently under progress. 
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